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Abstract. In this paper the relation between Pommaret and Janet bases of poly- 
nomial ideals is studied. It is proved that if an ideal has a finite Pommaret basis then 
the latter is a minimal Janet basis. An improved version of the related algorithm 
for computation of Janet bases, initially designed by Zharkov, is described. For an 
ideal with a finite Pommaret basis, the algorithm computes this basis. Otherwise, 
the algorithm computes a Janet basis which need not be minimal. The obtained 
results are generalized to linear differential ideals. 



1 Introduction 

Pommaret and Janet bases may be cited as typical representatives of invo- 
lutive bases |^J^,^,Q. Involutive bases are Grobner bases, though, generally 
redundant, and involutive methods provide an alternative approach to com- 
putation of Grobner bases. In so doing, polynomial Pommaret bases which 
were first introduced in ||^ have become a research subject in commutative 
algebra. They can be considered as generalized left Grobner bases in the 
commutative ring with respect to non-commutative grading Pommaret 
bases of homogeneous ideals in generic position coincide with their reduced 
Grobner bases Q . The use of these bases makes more accessible the struc- 
tural information of zero-dimensional ideals H]. Pommaret bases provide an 
algorithmic tool for determining combinatorial decompositions of polynomial 
modules and for computations in the syzygy modules . 

Linear differential Pommaret bases form the main tool in formal the- 
ory of linear partial differential equations ]To| , pd] | whereas linear differential 
Janet bases form an algorithmic tool in Lie symmetry analysis of nonlinear 
differential equations ||l^,|l^. Unlike reduced Grobner bases, Pommaret and 
Janet bases along with any other involutive bases lead to explicit formulae 
for Hilbert functions and Hilbert polynomials - 
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Basic properties of Pommaret and Janet bases are determined by the un- 
derlying involutive divisions . Non-noetherity of Pommaret division [Q| 
is responsible for non-existence of finite Pommaret bases for some polyno- 
mial (linear differential) ideals of positive (differential) dimension. On the 
other hand, any polynomial ideal as well as any linear differential ideal has 
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a finite Janet basis due to the noetherity of Janet division. Tlie two divi- 
sions differ greatly by their definition: unlike Janet divisibility, Pommaret 
divisibility do not depend on the leading terms of generators. Given an ideal 
and an admissible monomial ordering, or ranking in the differential case, its 
monic Pommaret basis is unique whereas there are infinitely many different 
monic Janet bases and among them only the minimal Janet basis is uniquely 
defined 

However, in spite of the above differences, Pommaret and Janet bases 
are closely related, and Zharkov was the first to observe this fact in the last 
paper of his life. He argued that if a polynomial ideal has a finite Pom- 
maret basis and a Janet basis which is Pommaret autoreduced, then they 
have identical monic forms (c.f. |^). Zharkov put also forward an algorithm 
for construction of Janet bases by sequential treatment of Janet nonmulti- 
plicative prolongations followed by Pommaret autoreduction. 

The goal of this paper is to study the relation between polynomial Pom- 
maret and Janet bases in more details and to improve the Zharkov algorithm. 
Our analysis is based on the properties of Janet and Pommaret divisions and 
involutivc algorithms studied in ||l|,p|,^Hl^ . 

This paper is organized as follows. The next section sketches some defini- 
tions, notations and conventions which are used in the sequel. Section 3 deals 
with analysis of the relationships between polynomial Pommaret and Janet 
bases. In particular, we prove that if an ideal has a finite Pommaret basis, then 
it is a minimal Janet basis. We describe here an algorithm of the combined 
Pommaret and Janet autoreduction which, given a Janet basis, converts it 
into another Janet basis. Since a minimal Janet basis is both Janet and Pom- 
maret autoreduced, the existence of a finite Pommaret basis is equivalent to 
Pommaret- Janet autoreducibility of any non-minimal Janet basis into a mini- 
mal one. In Section 4 we describe an algorithm for computation of polynomial 
Janet bases which is an improved version of the Zharkov algorithm jl^ . One 
of the improvements is the use of Pommaret- Janet autoreduction rather then 
the pure Pommaret autoreduction. Another improvement is incorporation 
of the involutive analogue [Q of Buchberger's chain criterion jl^. Section 5 
contains generalization of the results of Sections 3 and 4 to linear differen- 
tial ideals. The generalization is based on paper |Q where general involutive 
methods and algorithms of papers are extended from commutative 

to differential algebra. 

2 Basic Definitions and Notations 

Let N be the set of nonnegative integers, and M — {xf^ ■ ■ ■ xf^ \ dt £ N} be 
a set of monomials in the polynomial ring M = K[xi, . . . , Xn] over a field K 
of characteristic zero. 

By deg{u) and degi{u) we denote the total degree of u G M and the 
degree of variable Xi in u, respectively. If monomial u divides monomial v we 
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shall write u\v. Throughout the paper we restrict ourselves with admissible 
monomial orderings JTsf >- which are compatible with 

Xi y X2 y ■ ■ ■ y Xn ■ (1) 

The leading monomial of the polynomial f E M. with respect to >- will be 
denoted by lm{f). If F C M is a polynomial set, then by lm{F) we denote the 
leading monomial set for F, and Id{F) will denote the ideal in R generated by 
F. The initial ideal of an ideal / e M with respect to the monomial ordering 
y will be denoted by iny{I). The support of a polynomial /, that is, the 
set of monomials occurring in / with nonzero coefficients will be denoted by 
supp{f). For the least common multiple of two monomials w, u G M we shall 
use the conventional notation lcm(u,v). 

Definition 1. |p^ , |l4[ For a monomial u = x'l^ ■ ■ ■ x'fl' with dk > the vari- 
ables Xj,j > k are considered as Pommaret multiplicative or P— multiplicative 
and the other variables as P—nonmultiplicative. For u = 1 all the variables 
are P— multiplicative. 

Definition 2. iQ Let t/ C M be a finite monomial set. For each 1 < i <n 
divide U into groups labeled by non- negative integers di , . . . , : 

[di, . . . ,di] ~ { u E U \ dj ^ degj{u), 1 < j < i }■ 

A variable Xi is called Janet multiplicative or J— multiplicative for u G [/ if 
i = 1 and degi{u) = max{degi(u) | v G [/}, or if i > 1, u G [di, . . . , c?i_i] and 
degi{u) = inax{degi{v) \ v G [di, . . . , di-i]}. If a variable is not J— multipli- 
cative for M G [/, it is called J —nonmultiplicative for u. 

For a polynomial / G M the Pommaret separation of variables into mul- 
tiplicative and nonmultiplicative is done in accordance with Definition |^ 
where u = lm{f). Analogously, for an element / G -F in a finite polyno- 
mial set F C M the Janet multiplicative and nonmultiplicative variables are 
determined by Definition § with u — lm{f) E U = lm{F). 

We denote by Mp(/), NMp{f) and by M.,{f,F), NMj{f,F), respec- 
tively, the sets of F— multiplicative, F— nonmultiplicative and J— multipli- 
cative, J— nonmultiplicative variables for /. A set of monomials in F— mul- 
tiplicative variables for u and J— multiplicative variables for u £ U will be 
denoted by F(u) and J{u, [/), respectively. 

Remark 3. The monomial sets P{u) and J(u, U) for any u, U such as u G C/ 
satisfy the following axioms 

(a) If G F(u, U) and v\w, then v G L{u, U). 

(b) liu.v eU and uL(u, U) fl vL{v, U) ^ 0, then 
u G U) 01 V £ uL{u, U). 

(c) liveU and v G uL{u, U), then U) C F(w, U). 

(d) liV QU, then L{u, U) C L{u, V) for aU ueV. 
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if one takes either P{u) or J{u,U) as L{u,U). The axioms characterize an 
involutive monomial division, a concept invented in [Q. Every monomial set 
L(u, U) satisfying the axioms generates an appropriate separation of variables 
into (L— )multiplicative and (L— )nonmultiplicative. As this takes place, an 
element u G ?7 is an L— divisor of a monomial w G M if w/u G L{u, U). In 
this case w is L— multiple of u. Using the axioms, a number of new divisions 
was constructed [^,0 which may be also used for algorithmic computation 
of involutive bases. 

All the next definitions in this section are those in ||^,|^ specified to Pommaret 
and Janet divisions. 

Definition 4. Given a finite monomial set C/, its cone C(U), P—cone Cp{U) 
and J— cone Cj(U) are the following monomial sets 

C{U) =UueuuM, Cp{U)^\JueuuP{u), Cj{U)^\JueuuJ{u,U). 



Definition 5. Given an admissible ordering a polynomial set F C M 
is called Pommaret autoreduced or P — autoreduced if every f F has no 
terms P— multiple of an element in lm{F) \ lm{f). Similarly, a finite poly- 
nomial set F is Janet autoreduced or J — autoreduced if each term in every 
f G F has no J— divisors among lm{F) \ lm{f ). A finite set F will be called 
Pommaret- Janet autoreduced or P J — autoreduced if it is both P— autoreduced 
and J— autoreduced. 

Remark 6. From Definition ^ it follows that any finite set U of distinct mono- 
mials is J— autoreduced 

(Vu, v<gU) [u^v) [ uJ{u, U) n vJ{v, C/) = ] . 



Definition 7. Given an admissible ordering >- and a polynomial set P C M, 
a polynomial G M is said to be in P —normal form modulo F if every term 
in h has no P— divisors in lm{F). Similarly, if all the terms in h have no 
J— divisors among the leading terms of a finite polynomial set P, then h is 
in J— normal form modulo P. 

A general involutive normal form algorithm is described in and an invo- 
lutive normal form of a polynomial modulo any involutively autoreduced set 
is uniquely defined. We denote by NFp{f,F) and NFj{f,F), respectively, 
P— normal and J— normal form of polynomial / modulo P. 

Pommaret or Janet autoreduction of a finite polynomial set P may be 
performed similar to the conventional autoreduction jl^Jl^. If H is ob- 
tained from P by the conventional, or J— autoreduction we shall write H = 
Autoreduce{F) or H = Autor educe, j{F), respectively. 
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Definition 8. A P—autor educed set F is called a Pommaret basis (P— basis) 
of the ideal Id{F) generated by F if 

(V/ e F) (Vx e NMpU)) [NFp{f ■x,F)^Q]. (2) 

Similarly, a J— autoreduced set F is called a Janet basis (J— basis) of Id{F) 
if 

(V/ G F) [yx G NMjif, F)) [ 7Vi^j(/ ■x,F)^0]. (3) 



In accordance with Definition [7| the nonmultiplicative prolongation / • x with 
the vanishing P— or J— normal form modulo polynomial set F — {/i, • • • , /,„} 
can be rewritten as 

m 

f ■ X = fi hi, hi G M, lm{f) ■ X y lm{fi hi) (1 < J < m) 

i=l 

where supp{hi) C P{fi) or supp{hi) C J{fi,F), respectively, for every poly- 
nomial product fi hi . 

Let Gp and Gj be Pommaret and Janet bases of an ideal /, respectively. 
Then, from Definition || it follows that 

he I iff NFp{h,Gp) = and he I iff NFj{h,G.j) = . (4) 

This implies, in particular, the equalities 

Gp{lm{Gp)) = C{lm{Gp)) , Cj{lm{Gj)) = C{lm{Gj)) . (5) 

It is immediate from §) that lm{Gp) and /m(Gj) are P— and J— bases of 
the initial ideal iny{I). 

Corollary 9. If for a P — autoreduced set Gp the equality of its cone and 
P—cone holds and lm{Gp) is a basis of the initial ideal in(Id{Gp))}, then 
Gp is a P~ basis of Id{Gp). Analogous statement holds for a J— basis. 

Whereas monic Pommaret bases much like to reduced Grobner bases are 
unique, every ideal, by property (^, has infinitely many monic Janet bases. 
Among them there is the unique J— basis defined as follows. 

Definition 10. A Janet basis G of ideal Id{G) is called minimal if for any 
other J— basis F of the ideal the inclusion lm{G) C lm{F) holds. 

Remark 11. Every zero-dimensional polynomial ideal has a finite Pommaret 
basis, and for a positive dimensional ideal the existence of finite Pommaret 
basis can be always achieved by means of an appropriate linear transforma- 
tion of variables . 
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3 Relation Between Polynomial Pommaret and Janet 
Bases 

Given a finite monomial set J7, Definitions |^ and ^ generally give different 
separations of variables for elements in U. 

Example 12. U = {xiX2, X2X3, x'^}. 



monomial 


Pommaret 


Janet 


Mp 


NMp 


Mj 


NMj 


X1X2 


X2,X3 


Xl 


Xl,X2, X3 




X2X3 


X3 


Xl,X2 


X2,X3 


Xl 


x^ 


X3 


Xi,X2 


X3 


Xi,X2 



Here is, however, an important relation between Pommaret and Janet sepa- 
rations: 

Proposition 13. |l^(see also Q). If a finite monomial set U is P—autore- 
duced, then for any u £ U the following inclusions hold 

Mp{u, U) C Mj{u, U), NMj{u, U) C NMp{u, U) . 



For U in Example |l^ the minimal Janet basis Uj and the Pommaret basis 
Up the monomial ideal Id{U) are 

Uj^U\J{xixl}, Up^Uj\JT=2WiX2}yJf=2Wixl}\JT=2{4x3}. (6) 

Note, that Up is infinite and Uj C Up. Below we show that the inclusion 
Gj C Gp holds for any minimal Janet basis Gj and Pommaret basis Gp if 
both of them are monic and generate the same ideal. Furthermore, the proper 
inclusion Gj C Gp holds iff P is infinite. 

The following algorithm, given a finite polynomial set _F e M and an 
admissible ordering y, performs PJ— autoreduction of F and outputs a PJ- 
autoreduced set H. In this case we shall write H = AutoreducepjF. 

Since the involutive P— and J— reductions which are performed in the 
course of the algorithm form subsets of the conventional reductions , the 
algorithm terminates. Furthermore, the while-loop generates the P— autore- 
duced monomial set lm{H). In accordance with Remark^ the Janet autore- 
duction in line 12 docs not affect the leading terms, and, hence, the output 
polynomial set is both Pommaret and Janet autoreduced. 
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Input: e M, a finite set; ^, an admissible ordering 
Output: H = Autoreduce p J (F) 



begin 

G:=%- H:=F 1 

repeat 2 

H ■=H 3 

while exist h ^ H such that 4 

lm{h) (E Cp{lm{H\{h})) do 5 

H:=H\ {h}; G := G U {h} 6 

for each g G G do 7 

G\{g}; f:^NFj{g,H) 8 

if / ^ then 9 

H:^HU{f} 10 

until if = ff 11 

ii := Autoreduce J {H) 12 
end 



Proposition 14. // algorithm Pommaret-JanetAutoreduction takes a 
Janet basis F as an input its output H is also a Janet basis of the same 
ideal, and H Q F. 

Proof. Let F be a Janet basis of the ideal Id{F) and H be a polynomial set 
which computed by the algorithm. Apparently, Id{H) = Id{F). Consider U = 
lm{F). If U is P— autoreduccd, then H initiated as F in line 1 does not change 
in the process of the algorithm, and, hence, H = F = Autor educe pj{F). 
Otherwise, consider the output polynomial sets H . Denote lm{H) by V . 
Then, by construction, V C U. 

Consider a monomial t G Cj{U) and show that t G Cj{V). Let m G C/ be a 
J— divisor oit, that is, t G uJ{u, U), and u G be such that v\u. If u — v, by 
property (d) of Janet division in Remark we are done. Let now u G U \ V . 
The while-loop provides that u G vP{v). li v = xf^ '''^t' with dk > 
(1 < fc < n), then from Definition |l| it follows that u = x^^ ■ ■ ■ x'^^^'^'^ ■ ■ ■ x^" . 

We have to prove that any variable Xi G J{u, U) {1 < i < n) satisfy 
Xi G J{v,V). Consider two alternative cases: i < k and i > k. In the first 
case, by Definition ^, both u,v belong to the same group [di, ■ ■ ■ , di^i] of 
monomials in U. It follows that Xi G J{v, V). In the second case, by Definition 
0, Xi G P{v) and Proposition 13 we find again that Xi is J— multiplicative for 



V as an element in V. 

Therefore, F is a Janet monomial basis of Id{U). Thus, by Corollary ^, 
is a J— basis. In so doing, every J—normal form computed in line 8 of the 
algorithm apparently vanishes. This implies the inclusion H C F. □ 



Corollary 15. A minimal Janet basis is Pommaret autoreduced. 
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Proof. Let G be a minimal Janet basis. From Proposition and Definition 10 



it follows that lm{G) is P— autoreduced. Thus, G is P— autoreduced. □ 

It is clear that, given a Janet basis, its PJ— autoreduction yields, generally, 
more compact basis than the pure P— autoreduction. 

Example 16. Consider polynomial set 

F — {xyzt — xz, xyz + z^, xzt + x^ , xy + z, zt + x} 

which is a Janet basis of the ideal Id{F) with respect to the degree-reverse- 
lexicographical ordering >- such that x ^ y >~ z >- t. Given F and >- as 
input, the algorithm Pommaret-JanetAutoreduction(P, >-) outputs the 
minimal J— basis 

Gi = {xzt + x"^ ,xy + z, zt + x} . 

If one uses the Pommaret normal form computation in line 7 instead of that 
of Janet, then the algorithm leads to the output 

G2 = {xzt + x^ , z'^t + xz, xy + z, zt + x} = Gi U {z'^t + xz} . 

The following theorem is the main theoretical result of the present paper and 
forms a basis of an algorithm for construction of Janet and Pommaret bases 
described in the next section. 

Theorem 17. Given an ideal I <Z R and an admissible monomial ordering 
>- compatible with the following are equivalent: 

(i) I has a finite Pommaret basis. 

(ii) A minimal Janet basis of I is its Pommaret basis. 

(Hi) If F is a Janet basis of I, then G — Autoreducej p{F) is 
its Pommaret basis. 



Proof, (i) => (iii): Suppose G = {gi, . . . ,gm,} which, by Proposition is 
also a J— basis of /, is not its P— basis. Our goal is to prove that a Pommaret 
basis of / is an infinite polynomial set. From Corollary || it follows that 

{3g e G) e NMp{g)) [ lm{g) -x^^Gp (lm{G)) ] . (7) 

Among nonmultiplicative prolongations g ■ Xk satisfying choose one with 
the lowest Imia) ■ Xk with respect to the pure lexicographical ordering )^Lex 
generated by (|l|). If there are several such prolongations choose that with the 
lexicographically lowest Xk, that is, with the lexicographically highest g. This 
choice is unique since G is J— autoreduced. 

We claim that Xk £ Mj{g,G). Assume for a contradiction that Xk 6 
NMj{g,G). Then from Janet involutivity conditions (||) we obtain 



(3/ e G) [ lm{g) ■ xu - lm{f) ■w\w^J {lm{f),lm{G)) ] . (8) 



On the Relation Between Pommaret and Janet Bases 



9 



In accordance with condition d?) w contains P— nonmultiplicative variables 
for / and from (||) it follows |^ that / ^Lex g. If w — Xj, then xj >- 
Xk, and both lm{f) ■ Xj and lm{g) belong to the same monomial group 
[ degi{g), . . . , degk-i{g) ] appearing in Definition ^. Hence, Xj € NMj{f, G) 
in contradiction to (^). Therefore, deg{w) > 2 and (H) can be rewritten as 

lm{g) ■ Xk = {lm{f) ■ v) ■ Xj 

with w — V ■ Xj. Show that lm{f) ■ v ^ P {lm{f)). If u G Cp {lm{G)) we are 
done. Otherwise there is Xjn\v such that Xm G NMp {lm{f)). Since lm{f) ■ 
Xm -<Lex " a^fc, our choicc of g and Xk implies the existence gi G G such 

that / • Xm = gi ■ vi and v/xm G P {lm{gi)). If vi = v/xm S P {lm{gi)) 
we are done. Otherwise we select again a P— nonmultiplicative variable for 
gi occurring in vi and rewrite the corresponding prolongation in terms of 
its P— divisor lm{g2) G lm{G). Continuity of Pommaret division 0] provides 
termination of the rewriting process with an element in g E G \ {g} such 
that lm{g) is a P— divisor of lm{f) ■ v. Because G is P— autoreduced, by 
Proposition |l^ lm{g) is also a J— divisor of lm{f) ■ v. By this means there 
are two different Janet divisors lm{f) and lm{g) of the same monomial that 
contradicts Remark || and proves the claim. 

Let now iJ be a P— basis of Id{G). Denote lm{g) ■ Xk hy u and show 
that u e lm{H). Suppose there is an element h e H \G, such that u = 
lm{h) ■ V with v G P{lm{h)). Then h -<Lex u and there is g G G satisfying 
lm{h) = lm{q) -w where w ^ P {lm{q)). Thus, there is a P— nonmultiplicative 
prolongation q ■ Xj with Xj\w such that lm{q) ■ Xj ~<Lex u and lm{q) ■ xj ^ 
Cp {lm{G)) that contradicts the above choice of g and Xk- 

Now consider monomial set U = lm{G) U {u] C lm{H). By Defini- 
tion H, Xk e NMp{u) and u ■ Xk is obviously the lexicographically lowest 
P— nonmultiplicative prolongation of elements in U . It is easy to see that 
u- Xk ^ Gp{U). Indeed, since Xk G Mj{g, G), it follows that u ■ Xk ^ U, and 
a P— divisor of u ■ Xk would also P— divide lm{g) that is impossible as G is 
P— autoreduced. Thus, we find that u ■ Xk G lm{H). By sequential repetition 
of this reasoning for it • a;^ (i = 2, 3 . . .) we deduce that every such monomial 
is an element in lm{H), and, therefore, H is infinite. 

(Hi) =^ (m): If G is a minimal Janet basis of Id{G), then Corollary 
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implies G = Autor educe jp{G), and the above arguments show that either G 
is also a P— basis or the latter is infinite. 

(a) (i): This implication is obvious. □ 



Corollary 18. Let Gj be a monic minimal Janet basis and Gp be a monic 
Pommaret basis for the same polynomial ideal and monomial ordering. Then 
Gj ^Gp and Gj C Gp iff Gp is infinite. 



Proof. This follows from the above proof of Theorem |l 



□ 
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Corollary 19. // a polynomial ideal is in generic position, then its minimal 
Janet basis is also a Pommaret basis. If such an ideal is homogeneous, then 
these bases are also reduced Grobner bases of the ideal. 

Proof. I has a finite Pommaret basis [HQ]. If / is homogeneous, then its 
Pommaret basis is the reduced Grobner basis 0. □ 

4 Algorithm for Construction of Janet Bases 

In this section we present an algorithm for constructing J— bases of poly- 
nomial ideals which is based on Theorem |l^ and will be called JanetBasis. 
Whenever the ideal generated by an input polynomial set has a finite P— basis 
for a given admissible ordering, the algorithm outputs just this basis which is 
also a minimal J— basis. Otherwise, the J— basis computed by the algorithm 
is not necessarily minimal as we demonstrate below by the explicit example. 



Algorithm: JanetBasis(F, )~) 



Input: P e M, a finite set; an admissible ordering 






Output: G, an involutive basis of the ideal Id{F) 






begin 






G := Autoreduce{F) 




1 






2 


for each g & G do 




3 


T:=TU{{g,lmig),9)} 




4 


while exist {g, u, D) G T and x e NMj{g, G)\D 


do 


5 


choose such {g, u,D),x with the lowest lm{g) 


• X w.r.t. >- 


6 


T:=T\{{g,u,D)]\j{{g,u,DVj{x})} 




7 


if Criterion{g ■ x, u, T) is false then 
h:^NF,j{g-x,G) 




8 




9 


ii h^O then 




10 


if lm(h) — lm{g ■ x) then 




11 


T:^TU{(/i,u,0)} 




12 






13 


T:=TU{{h, lm{h),(i))} 




14 


G := AutoreducepjiG U {h}) 




15 






16 


T 




17 


for each q ^ G do 




18 


if exist (f,u, D) e Q s.t. lm(f) = lm(q) then 


19 


choose g € G s.t. u G lm{g)J (lm{g), Im 


(G)) 


20 


T := T U {(.g, /m(g), D n NMj{g, G))} 




21 


gIsg 




22 


T:=TU{(g,M5),0)} 




23 


end 







Criterion{g,u,T) is true provided that if there is {f,v,D) £ T such that 
lcm{u,v) -< lm{g) and lm{g) £ lm{f)J {lm{f),lm{G)). 
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The structure of algorithm JanetBasis is very close to that of the algo- 
rithm InvolutiveBasis devised in for construction of involutive bases for 
arbitrary constructive involutive divisions, and, hence, applicable to Janet 
division. The main difference between the algorithms is the form of their in- 
termediate autoreduction. Whereas the previous algorithm, when specified 
for Janet division, uses the pure Janet autoreduction which do not affect 
the leading terms, the below one uses the above algorithm Pommaret- 
JanetAutoreduction. By this reason, the algorithm InvolutiveBasis, un- 
like the below one, almost never outputs a minimal Janet basis or a Pommaret 
basis if the latter is finite. In paper we designed the algorithm Minimal- 
InvolutiveBasis which always outputs a minimal involutive basis whenever 
the latter exists. As we now see from Theorem 0, this algorithm in the case 
of Janet division outputs also a Pommaret basis if it is finite. 

Besides, in the algorithm InvolutiveBasis the involutive autoreduction 
is performed whenever nonzero normal form is obtained in the while-loop. 
In the algorithm JanetBasis the subalgorithm Pommaret-JanetAutore- 
duction is caused by a nonzero J— normal form h in line 15 only if the leading 
term of the related prolongation is J— reducible. Otherwise, since G is always 
P— autoreduced before its enlargement with h, lm{h) cannot P— divide, in 
accordance with Definition 0, any other element in lm{G). 

Note that we indicated the intersection in line 21 to emphasize that ele- 
ments in D must be nonmultiplicative variables for the corresponding poly- 
nomial that is always understood in algorithms of papers ( Jl|,pPJT^ ) . 

Noetherity of Janet division provides termination of the algorithm Janet- 
Basis 10]. To show this consider the intermediate bases Go = Autoreduce{F) 
and Gi (i — 1,2,...) generated after the i — th iteration of the while- loop. 
It is clear that 

Id (/m(Go)) C Id (lm{Gi)) C Id (Zm(G2)) C • ■ • (9) 

and this chain is stabilized after finitely many steps. Namely, the stabilization 
starts when the intermediate polynomial set G becomes a (non-necessarily 
reduced) Grobner basis of Id{F). By partial involutivity of G Q, the proper 
inclusion in chain (^) holds only when G is enlarged by h in line 15. In between 
of such proper inclusions and after the chain stabilization lm{G) is completed 
with lm{h) — lm[g ■ x) as stands in line 11, and this completion cannot be 
infinite by noetherity of Janet division ||l[. Once algorithm terminates, it 
produces, by Proposition a PJ — autoreduced Janet basis of F because 
the involutivity conditions (H) are satisfied as is checked in lines 6 and 9 where 
h = 0. In so doing, correctness of the criterion which is verified in line 8 is 
proved exactly as it done in [|l| for the algorithm InvolutiveBasis. 

Remark 20. In line 6 of the algorithm JanetBasis one can use any admissible 
ordering for selection of the current J— nonmultiplicative prolongation g ■ x 
to be treated in the following lines. This selection ordering may not only be 
different from the main ordering >~ but also may vary at every step when the 



12 Vladimir P.Gerdt 



selection is done. Correctness of this arbitrariness in the choice of selection 
ordering follows from the correctness of this arbitrariness for the monomial 
completion procedure JT6| . 

As mentioned above, the algorithm JanetBasis for an ideal of positive di- 
mension may not output its minimal J— basis. We demonstrate this by the 
following example. 

Example 21. Let be a set {x^y — z,xy^ — y} generating one-dimensional 
ideal, and >- be the degree-reverse-lexicographical ordering with x )~ y ^ z. 
Then the algorithm JanetBasis (i^, >-) outputs the following Janet basis 

( 2 23 22 2-1 

G = \x y — z,x z — z,xy — yz,xz — z ,y z — y,y z — z) 

whereas the minimal Janet basis coinciding with the reduced Grobner basis 
is 

{xy - y z,xz - z'^.y'^z - y,y z"^ - z) . 



Remark 22. The algorithm JanetBasis is an improved version of the algo- 
rithm designed by Zharkov ||l5|] . The first improvement is the use of the mixed 
Pommaret- Janet autoreduction instead of the pure Pommaret autoreduction 
as Zharkov proposed. Let d and G2 be output Janet bases computed with 



the use of PJ— and P— autoreduction, respectively. Then, Proposition [13 
implies Gi C G2 and below we give an example when the proper inclusion 
holds. The second improvement is the criterion used in line 7. This criterion 
is an involutive analogue Q of the Buchberger's chain criterion ||l^ and is 
superior to the criterion used in p5[. 



Example 23. The following polynomial set generates three-dimensional 
ideal 



4x5(x6 + 2ai — 8xi)(a2 - 03) - X2X3X4 + X2 + X4 
4:X5{x6 + 2ai — 8a;2)(a2 - 03) - X1X3X4 + a;i + 2:3 
4x5(a;6 + 2ai - 8a;3)(a2 - 03) - xiX2Xi + X2 + x^ 
Ax^ixe + 2ai - 8a;4)(a2 - 03) - X1X2X3 + xi + X3 



where 

ai = xi + X2 + X3 + X4, a2 = X1X2X3X4, 03 = X1X2 + X2X3 + X3X4 + X4X1 



For the polynomial set and the degree-reverse-lexicographical ordering 
compatible with (|^) the algorithm Janet basis outputs set Gi with 71 poly- 
nomials. The pure Pommaret autoreduction in line 15 generates set G2 D Gi 
of 75 elements. Note that a minimal Janet basis contains 49 polynomials 
whereas the reduced Grobner basis contains 44 polynomials. 
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5 Linear Differential Bases 

Let K be a zero characteristic differential field with a finite number of mu- 
tually commuting derivation operators d/dxi, . . . , d/dxn- Consider the dif- 
ferential polynomial ring BM = IC{yi, . . . ,ym} with the set of differential 
indeterminates {j/i, . . . , Hm}- Elements in DR are differential polynomials in 
{Ui, . . . jUm}- An ideal in DM generated by linear differential polynomials is 
called linear differential ideal . 

In Ql , by exploiting the well-known algorithmic similarities between poly- 
nomial and linear differential systems and the association between monomials 
and the derivatives, we extended the general involutive methods and algo- 
rithms designed in Jl|,p|jT^] for polynomial ideals to linear differential ideals. 
The statements and algorithms of Sect. 3 and 4 admit similar extension. As 
this takes place, all the above statements have proven for the polynomial case 
are proved by parallel arguments for the differential case. In the following 
table we give a short correspondence between these two cases. In particu- 
lar, this correspondence allows one to rewrite the algorithms Pommaret- 
JanetAutoreduction and JanetBasis for linear differential bases. 



Commutative Algebra 


Differential Algebra 


R= K[xi,...,Xn] 


DM = IfC{j/i,...,y„} 






F,G,H cR 


F,G,H c DM 


X = Xi ■ x„ 




g ■ X 


dxg 


monomial ordering 


ranking 


lm{g) 


leading derivative ld{g) 


lm{G) 


ld{G) = UgeclWCg)} 



6 Conclusion 

As we have seen finite Pommaret bases of polynomial and linear differential 
ideals are minimal Janet bases. The above proof of Theorem [T7| shows that, 
given a P— autoreduced Janet basis G, it is easy to verify the existence of 
a finite Pommaret basis. One suffices to check the condition (^. Another 
check which may be even easier in practice is to look at the Pommaret and 
Janet separation of variables for elements in G. As shown in Q, the exis- 
tence of a finite Pommaret basis implies coincidence of both separations for 
every element in G. Moreover, the leading monomial structure of an infinite 
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Pommaret basis can be read off the structure of lm[G) (c.f. (g) for Example 
©■ 

Therefore, minimal Janet bases can be used in commutative and differ- 
ential algebra as well as a Pommaret bases with the advantage of finiteness. 
For example, in the formal theory of differential equations |^,|l^,|ll| infinity 
of a Pommaret basis signals on (5-singularity of the coordinate system chosen, 
and the condition (^) gives the same signal for Janet bases. 

On the face of it, Pommaret division looks like more attractive than 
Janet one since its separation is easier to compute than the Janet separa- 
tion. Besides, Pommaret division, unlike that of Janet, is globally defined 
Hence, after enlargement of the intermediate polynomial set with an 
irreducible P— nonmultiplicative prolongation there are no needs to recom- 
pute the Pommaret separation for other elements in the set. However, the 
careful implementation of both divisions do not reveal any notable advan- 
tage of Pommaret division over Janet division in construction of polynomial 
bases. This rather surprising fact was firstly observed by Zharkov One 
of the explanations of this experimental phenomenon is given by Proposition 

one must generally treat more P— nonmultiplicative prolongations than 
J— nonmultiplicative ones. In so doing, the search for an involutivc divisor in 
the process of involutivc reduction can be done similarly for both divisions 
as we show in our forthcoming paper pT[ | . 

The algorithm JanetBasis presented above is now under implementation 
in C in parallel with the algorithm MinimallnvolutiveBasis |^ specified for 
Janet division. Our first experimentation with the codes shows that some- 
times the former algorithm runs faster than the latter one and needs less 
computer memory. For example, modular computation of the degree-reverse- 
lexicographic Janet basis for the Cyclic 7 example ||2^ is about twice faster 
with the algorithm JanetBasis than with the algorithm Minimallnvolu- 
tiveBasis. Currently, the timings for computation modulo 31013 are about 
5 and 10 minutes, respectively, on a Pentium Pro 333 Mhz computer. 
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